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Synopsis : The transient network model of Green-Tobolsky [1946], Yamamoto [1956] and Tanaka- 

Edwards [1992] is formulated within the frame of thermodynamics of irreversible processes, using as a 
fundamental quantity the chemical potential associated to the connection of strands to the network and treating 
these connections as chemical-like reactions. All thermodynamic quantities are thus naturally defined in and 
out of equilibrium. Constitutive equations are derived, giving the stress and the heat production as functions of 
the thermomechanical history. The Clausius-Duhem inequality, stating that the source of entropy is non- 
negative, is shown to hold for any thermomechanical history, ensuring the thermodynamic consistency of our 
model. The presented model includes the Green-Tobolsky model, whereas those of Yamamoto and Tanaka- 
Edwards fit within ours on the condition that their free parameters obey a detailed balance condition stemming 
form Boltzmann equilibrium statistics. 
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Introduction 

The transient network models first aimed at a molecular description of the viscoelasticity of polymers in 
the entangled regime, either in the melt state or in concentrated solutions. These systems exhibit a rubberlike 
elasticity when subjected to fast deformations but, in contrast with rubber, their stress relaxes when the 
deformation is kept constant. The pioneering work of Green and Tobolsky (1946) started the subject by 
suggesting that the laws of rubber elasticity were valid in non-crosslinked polymeric systems, the rubber 
crosslinks being replaced by the entanglements between chains. An important characteristics of entanglements 
is that Brownian motion enables the chains them to disentangle and reentangle with other chains, in contrast 
with the permanent character of the rubber chemical crosslinks. This picture of entangled systems has been 
since modified by the reptation theory, but a more definite embodiment of transient networks is provided by 
associating polymers, which carry sticking sites able to form transient bounds with one another. There exists a 
great variety of binding mechanism, such as solvent complexation, local crystallisation or hydrogen bonding, to 
name just a few. Recent interest has been focused on water solutions of end-capped associating thickeners, 
obtained by grafting hydrophobic paraffinic groups at the ends of a water-soluble polymer chain. In solution, 
these strands join their end groups in tiny micelles acting as the connection points of a three-dimensional 
network. Brownian agitation and/or mechanical action can force an end group out of a micelle and, after 
having diffused in water, this end group will eventually join another micelle. The elastic links between micelles 
are provided by the hydrophilic chains, which can be made monodisperse, making this system a well- 
characterised and simple transient network [Annable et al (1992)]. 

All transient networks share the following characteristics: They are made of elastic strands joined by 
temporary, reversible junctions. The creation or the disappearance of junctions causes the strands to connect to 
or to disconnect from the network. The junctions' temporary character causes the network to be continuously 
reworked, some connected strands disconnecting while free ones connect. When the sample is submitted to 
deformation, the disappearing network has been strained by the deformation accumulated during its lifetime, 
while the newly created one incorporates strands that carry less or no stress, having had time to relax before 
their connection. The stress relaxation thus directly originates in the transient character of the network, and the 
relaxation time is related to the junction lifetime. 

Such simple characteristics lend themselves to a molecular description of both the thermodynamics and 
the viscoelasticity of transient networks. The first attempt is the Green-Tobolsky model, which considers a 
network made of Gaussian strands which are created at a constant rate according to an equilibrium distribution, 
and have a constant probability per unit time to disappear, i.e. to disconnect from the network. As a 
consequence of these assumptions, Green and Tobolsky (1946) showed that the stress is related to the strain 
history according to a codeformational Maxwell equation, and that isothermal periodic closed-cycle 
deformations lead to dissipation of mechanical energy. The codeformational Maxwellian stress-strain relation 
however suffers from limitations: the simple shear steady-state viscosity shows no shear-thinning, the second 
normal stress difference vanishes, and the elongational viscosity diverges at a critical extension rate. 
Yamamoto (1956-1957-1958) addressed these difficulties by introducing a strand disconnection rate that 
increases when the strand is extended, and by relaxing the Gaussian strand assumption. He was able to work 
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out an expression for the instantaneous energy dissipation. This model better agrees with experiment when a 
moderate dependence of the disconnection rate versus the strand extension is assumed, but when pushed too 
far, it can predict thermodynamically impossible behaviours. Yamamoto himself pointed out in his 1957 paper 
that the viscosity can be made negative by choosing a disconnection rate proportional to n-th power of the 
strand extension, with n>5. The same difficulties arise in the Tanaka-Edwards model (1992a-d), which adds to 
the preceding one a balance between connected and disconnected strands, making the connection rate 
proportional to the number of disconnected strands. 

Our aim is to formulate the simplest complete rheothermodynamic model of transient networks including 
Yamamoto and Tanaka-Edwards' ideas. We call « Rheothermodynamics » the thermodynamics of systems 
submitted to deformation, and « complete » means that all thermodynamic quantities, including the stress and 
the rate of heat production as well as the thermodynamic potentials, can be related to the thermomechanical 
history of the system, i.e. the past and present temperature and deformation. The model is kept as simple as 
possible by minimising the number of assumptions, with no recourse to special mechanisms which would only 
operate out of equilibrium. Since most thermomechanical histories drive the system out of equilibrium, our 
thermodynamic description must be framed within the non-equilibrium thermodynamics of irreversible 
processes. The present paper comprises seven sections arranged in order of increasing specialisation. The first 
section formulates the general thermodynamic relations who do not depend on the precise form of the 
relaxation phenomena, using as fundamental quantities the strand distribution function and its 
thermodynamically conjugate quantity, the strand chemical potential. The second principle is formulated in 
terms of the Clausius-Duhem inequality. The second section specialises to ideal systems, corresponding to non- 
interacting strands, and the third one examines the constraints the Boltzmann equilibrium statistics places on 
the model. The fourth section introduces the kinetic processes, namely the connection/disconnection reactions 
and the heat conduction, allowing to formulate the evolution equation for the distribution function, thus 
completing the rheothermodynamic model this paper aims at. The model is then shown to be 
thermodynamically consistent for any given thermomechanical history. The special case of entropic systems, 
such as polymers, is treated in the fifth section. The sixth section is devoted to the limiting case of a connection 
rate independent of the number of connected strands. Eventually, the seventh section considers some restricting 
assumptions usually made in the transient network literature and examines how the formerly published models 
fit within ours, uncovering some thermodynamic inconsistencies. 
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1/ General 



The system we consider is made of strands, in solution or in the melt, able to form a network through mutual 
transient connections. When connected, strands join two different points of the network separated by a vector 
h we call the connector. The connection to the network as well as the disconnection are instantaneous events, 
and the connection/disconnection rate is supposed to be slow enough for the kinetic degrees of freedom to 
equilibrate within this time scale. The network is constrained to deform affmely to the macroscopic 
deformation, whereas the free strands are assumed to undergo a Brownian diffusion in the h space. The 
fundamental point of this paper is that the connection/disconnections can be treated as chemical-like reactions 
between free and connected strands, and the motion in the h space can be considered as well as a continuous 
chemical reaction in the sense of Prigogine and Mazur (1953). With such provisions, transient networks belong 
to the field of thermodynamics of irreversible processes in systems with internal degrees of freedom [Prigogine 
and Mazur (1953), de Groot and Mazur (1984) § X-6]. The thermodynamic functions are thereby defined in 
non-equilibrium situations exactly as in a system undergoing a chemical reaction. For the sake of simplicity, we 
consider only one kind of strands ; the generalisation to a distribution of strand species can be made along the 
line of Lodge (1968), by adding the contributions of the various species to the various thermodynamic 
potentials. In the same spirit of simplicity, the strands we consider have only two possible connection states : 
they can be either free or connected, therefore dangling strands are ignored. 

Consider a sample of volume V, containing N strands and N s solvent molecules. This sample is taken small 

enough for space-dependent quantities such as the stress, strain and temperature to be considered uniform 

within V. The strands are characterised by a connection state a taking one of the two possible values F, for free, 

and C, for connected to the network, and by a connector h with the dimension of a length, belonging to the 

three-dimensional space tangent to the physical space. The strand distribution function r ^ is such that 

T r d 3 h is the number of strands within volume V, in connection state a, and with vector h within the 
ah 

differential volume element d h = dh l dh 2 dh 3 . Y ? has the dimension [length] , the reciprocal of a volume in 
the h -space, r r is thus an intensive variable in the h -space and an extensive variable in the physical space. 
The number N of strands is the sum of the numbers N a of strands in connection state a 

N = Y,N a ; N a =\d 3 hr ai (1) 
a 

Here and in the following, the sum ^ runs over the states a=F and a=C, and the integrals with respect to h 
extend over the whole h space. We now restrict our attention to systems homogeneous enough for the diffusion 
in the physical space to be negligible. Since in addition the solvent molecules, as well as the strands, are neither 
created nor destroyed, their respective numbers are constants: 

dN s =0 and dN = dN F +dN c =0 (2) 

The thermodynamic state of the sample is determined by its volume V, its temperature T, the distribution 
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function r ^ and the number N s - Since N s is a constant, the Gibbs equation for the internal energy E is written 
[Prigogine and Mazur (1953), de Groot and Mazur (1984)] 



dE = TdS - PdV + 



a 



where \i ^ is the chemical potential associated to an (a h ) -strand, S is the entropy, P is the pressure and T is 
the absolute temperature, fi ^ reduces to the usual chemical potential per molecule whenever strands identify 
with molecules. In any case, Ur is the quantity thermodynamically conjugate to Y r . Eq.(3) is written for 

IX / / IX / / 

independent variations dV and dT a ^, this means that two kinds of external forces are separately acting on the 
sample : first, the pressure, conjugate to the volume, and second, the forces acting on the network, changing the 
connector h and the distribution of connected strands. Note that the volume is a thermodynamic variable 
but the shape is not [Landau and Lifshitz (1980)], in that a transformation which changes the shape but brings 
V, T and T a ^ back to their original values will leave all thermodynamic functions unchanged. 
The macroscopic deformation must now be considered. Let x be the position vector of a given material point at 
time t, and let x be its value at time f. The deformation gradient is the tensor relating the differentials dx and 
dx : 

F(tAj=§; ; FM=FV,ty x (4) 

Throughout this paper, Cartesian indices will be denoted i, j, k and I, and will be submitted to the summation 
convention. The velocity gradient K is related to F(t,t') considered as a function of the present time t, the 
reference time t' being kept constant, according to 



_ dv t _ dF(t,t') ik 

Because of the absence of diffusion in the physical space, the network must, on average, follow the macroscopic 
deformation. The affine deformation assumption we now introduce states that connected strands deform 
affmely to the macroscopic deformation, their connector h at time t being linked to their value K at time f by 
the same relations as the differentials dx and dx 

dh- dF::(t,t') 

h i =F u (f,t')h' i and — '- = — ±- (6) 

dt at 

provided that the strand stays connected during time interval (t, f ')• The affine deformation assumption amounts 
to linking the two sets of forces previously considered by a coherency constraint in the sense of Sekimoto 
(1991), ensuring that the solvent does not flow through the network. In short, we consider the limit of strong 
friction between network and solvent. 

Macroscopic deformations thus change the distribution of connected strands. In addition, the distribution of 
both free and connected strands varies under the action of the kinetic relaxation mechanism to be introduced 
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later in section 4, namely the Brownian motion and the connection/disconnection reactions. The variation 
dT a f t in an infinitesimal deformation will therefore be split into two terms 



dr ci= dAr ci +dRr ci 



(7) 



where d R T a ^ accounts for the still unspecified relaxation mechanisms. The affine variation of the connector h 
amounts to a flow in the h -space, so d A T a ^ will be written 



d A T crt * i \ dC^ 



= -4-^ch h J K ij)= -^r h i K ij ~ r ch 



dt dhi VCh J ' " dhi 



Kick (8) 



The first term of the right member stems from the affine variation of the connector, and the second term 
accounts for the variation of the volume element d h . 
The extra-stress tensor T is defined as 



X ij=-h\ d ^ch 



i 



h j +8 a Y ch 



(9) 



and, assuming that T c ^ vanishes fast enough as h goes to infinity, an integration by parts yields 

l f r fy C h 

T .. =- \d 3 hT r r— — hj (10) 
V y) Ch dh , J 

The contribution of d A to the Gibbs equation (3) can then be written as 

\d'hn cti d A Y^=V, ijKij dt (11) 

and the Gibbs equation then becomes 

dE = TdS+VG ij K i jdt + Y,\d 3 hn ali d R T afi (12) 
a 

for an infinitesimal change taking place within time interval dt, where the total stress tensor 

a i j=-P5 ij + z ij (13) 

is the quantity conjugate to VK^dt . The quantity -PdjjVKjjdt = -PdV is the work performed by the pressure, 
and since VT^K^dt is the work performed in affmely deforming the h vectors, the extra stress T is the stress 
associated with the variable h . The stress a thus represents the total effect of the mechanical forces conjugate 
to both the network deformation and the macroscopic deformation of the sample, locked with one another by 
Sekimoto's coherency constraint resulting in the affine deformation assumption. 

Expression (13) is derived under very general hypotheses that do not imply that T is symmetric. A non- 
symmetric T can result, for instance, from the torque an electric field imposes on strands having an anisotropic 
polarisability, contributing an orientation-dependent term to their chemical potential. A sufficient condition 

6 



for the symmetry of f will be given in the seventh section. 

Introducing the heat flux density q allows the first principle to be written as 

dE 

= -V div q + Va ti K u (14) 

dt y y 

where the quantity Vdtdiv q is the heat flowing out of volume V during time interval dt. Volume V is taken 
small enough for div q to be practically homogeneous within V. This also holds for div (q/T) and q-VT~ l , so 
using eq.(12) and (14) permits to split dS into a flux term and a source term accounting for irreversible 
processes: 

dS .... q d irr S 

— = -Vdiv— + (15) 

dt T dt 



where the entropy source reads 



= Vq^---V\d 3 hfl r — (16) 

dt T T q ^ ah dt 

The term div (q/T) is the divergence of the entropy conduction flux, and q-VT is the rate of entropy 
production per unit volume due to heat conduction [de Groot and Mazur (1984)]. 

The irreversible phenomena, namely the heat flux q , the Brownian diffusion flux of free strands in the h 
space and the connection/disconnection reaction rate, obey a set of constitutive equations which is left 
undefined up to the fourth section of this paper, where a specific model will be formulated. However, in order to 
comply with the second principle, the constitutive equations must satisfy the Clausius-Duhem inequality 

^/>0 (17) 
dt 

ensuring that the entropy source is non-negative, whatever the thermomechanical history. A stronger 
requirement, valid in the case of Brownian diffusion [Prigogine and Mazur (1953), de Groot and Mazur 
(1984)], states that the integrated term of eq.(16) must be non-negative rather than the integral as a whole, so 
the entropy generated within the volume element d h is non-negative. This condition is satisfied by the kinetic 
equations presented in the fourth section of this paper. 

A direct consequence of inequality (17) is that the isothermal steady flow viscosity T}(jc) = C y K"y /fyfy must be 
non negative, provided that the system is able to reach a stationary state characterised by steady values of the 
thermodynamic variables. The inequality ?7(fC )> follows from the stationarity conditions dE/dt = and 
dS/dt = 0, inserted into eq.(15) and (16) respectively. The volume is a thermodynamic variable, therefore it 
must altered by the flow, hence the condition = . 

The thermodynamic completeness of our formalism permits to derive a heat equation. The heat released by the 
sample can be obtained from eq.(15), by inserting the differential dS evaluated as a function of T, V, and T ^ : 
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JV,{T) 



dS(T,V,{T ari }) = -C v{r] dT + 



( ^ah A 



dr 



JV,{T} 



where the subscript {r) means that the partial derivative is performed at constant distribution r^. The first 
term of the right member involves the heat capacity per unit volume, at constant volume and distribution, 

TfdS^ 



-V{T} 



JV,{T} 



(19) 



The second term follows from a Maxwell relation involving the free energy A(T,V, {T^ }) = E-TS : 



d 2 A 

dvdr 



dp_ 



(20) 



and another Maxwell relation allows writing the integrand of the third term as 



SS 

ST r 
ah 



Jry 



SdA 
Kh^ 



^ah 



Jv 



dr 



(21) 



where SS I ST a ^ denotes the functional derivative of S with respect to T n ^ [Courant and Hilbert (1953)] 



ah 



Finally, splitting dT a ^ according to eq.(7)-(8) yields 



^(r,y,{r a ^ }) = ^c y{r} jr-^ {r} ..^-£j^ 

a 



d ^ah 



dr 



dY ah ^ 



JV,{T} 



where we introduce , the entropic stress tensor at constant distribution, such that 



dr 



d T S TfdP^ 



dt 



y^ 1 nat) 



dV 
dt 



(23) 



s {r} is then the logarithmic derivative of the stress with respect to the temperature at constant volume and 
distribution : 



where the tensor 



s {T}ij ~ T 



dr 



= -T 



5 V +t V 



(24) 



t - = T 



-{d^hT-h:^- 

V i Ch 1 dhi 



dr 



(25) 



is accordingly named the entropic extra stress tensor. It will be shown in section 5 of this paper that systems 
that exhibit entropic strand elasticity are such that t = x , i.e. their extra stress is entirely entropic. 
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The rate of heat release can be written, from the expression (15) and (22) of dS, as 

diva = -C V r n — + 5' r „...K- ;; +— Y \d?h 

1 1 dt f r iy 1 V ^ 

If the deformation rate is faster than the relaxation mechanisms, then the factor d R T a ^jdt can be neglected in 
eq.(22) and (26), which then reduce to divq = -T/V dS/dt = -Cy^dTjdt + s^ytCy . The entropic stress thus 
governs the transformation of mechanical energy into heat in fast deformations. However, it may happen that 
some of the relaxation mechanisms are fast, at least faster than the deformation rate, like the Brownian 
diffusion of free strands in the model presented in section 4. The preceding interpretation will then be shown to 
hold with modified coefficients of dT/dt and K" . 



dT 



ah 



d R T ~ 
u L ah 

dt 



(26) 
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21 Ideal systems 



The chemical potential of strands must now be explicited as a function of the thermodynamic variables. The 
simplest case is that of non-interacting strands constituting an ideal system, where the chemical potential takes 
the form [Prigogine and Mazur (1953), de Groot and Mazur (1984) § X-6] 



er 



ah 



ah 



ah 



N 



(27) 



where k B is the Boltzmann constant, and the constant 9 is a volume of the h -space, making the argument of 
the logarithm non-dimensional. The standard chemical potential £ ^ is the value the chemical potential takes 
when the concentration T^l N amounts to one strand per volume 6 of the h -space. £ ^ thus accounts for 
the interactions of the strand with its surrounding molecules, whilst kgTlnOT^ j N is the combinatoric 
entropy of the strand distribution. In the following, f r will be called simply the connector potential. Since the 
system is ideal, f r depends on the intensive variables T and P (note that at constant T and P, the volume 
varies with the strand distribution, the strand volume V a ^ = \^fi a ^ /dpjj, depending on a and h . is 
thus a function of T, Vand of the whole distribution {r}). 

The gradient of the chemical potential in the h -space splits into two forces [Jongschaap & al., (1997)] 



3lnr 



ah 



dh 



(28) 



The first term is the connector force 



f ~=- 

J ah 



dh 



(29) 



deriving from the connector potential, and defined with the customary positive sign. The second term is the so- 
called the Brownian force, coming from the combinatoric entropy. Inserting expression (28) into the definition 
(10) of the extra stress tensor yields the familiar expression 



^J d3 * r c*(w^-*« r *#) 



(30) 



The entropic extra stress t , defined in eq.(25), has the same form as x in eq.(10) with 9/* a ^ j~dh substituted 
for 



d ^ah 



dT 



= T 



dT 



dlnr r 
+ k»T =A 



(31) 



The first term is the entropic part of the connector force and the second term is just the Brownian force, which 
is entirely entropic. The entropic extra stress then reads 



t ij = y\^ T Ch 



d f C hi 
dT 



hj-k B dij 



Jv 



(32) 
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3/ Equilibrium 



Since Brownian diffusion and connection/disconnections allow the strands to reach every possible state (a, h) , 
the chemical potential at equilibrium cannot depend on a nor h : 



9T e 



» e JPJ) = »:\=L rh {P,T) + k B T\n- 



" ah 



ah ^ah x - - » ^ 

Inserting this relation into (10) gives the equilibrium extra stress tensor: 



= eq 

T =0 



(33) 



(34) 



: eq 

There is no such general relation for the entropic extra stress at equilibrium t , because equation (33), stating 
that /u a ^ is uniform in the (a,h) -space, does not imply that \^dfi a ^ ^oT j also is uniform. An example of 
this situation is given about eq.(87). 

The equilibrium distribution function obeys the Boltzmann statistics 



1-1= ^exp 
ah Q 



' ^eq £ a h 



knT 



(35) 



We define for further use the normalised equilibrium distribution (j> a ^ , such that ^d 3 h(/) a ^ = 1 

~^ah 



T e 1 
ah 



exp 



N, 



eq r 

a \d 3 h' 



exp 



^ah' 



(36) 



It is linked to the connector force by the relation 

dln< t>ah 



fah=- k B T - 



dh 



(37) 



Introducing the connection potential = C, c ^ - C, F ^ as the difference of connector potentials with the same 



connector, and using the equilibrium condition (33), we rewrite the ratio T eq ~ T eq - of equilibrium 

Fh I Ch 

distributions according to the mass action law 



Fh 
yeq 
Ch 



= exp 



(38) 



In the important case of an h -independent connection potential, the equilibrium distribution of connected 
strands and that of the free strands have the same dependence on h . 
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4/ Kinetics 



The formulation of kinetic equations for the dissipative processes d R T a ^ and q permits to derive explicit 
rheothermodynamic constitutive equations, giving the thermodynamic functions, the distribution function, the 
stress and the heat production as functions of the deformation and temperature history. These kinetic equations 
must result in a non-negative source of entropy, in order to comply with the second principle. 
We first formulate the kinetic equations. The connection and the disconnection of strands are assumed to take 
place with no change of h , fast enough to be considered instantaneous at the thermorheological timescale. This 
allows treating the connections and the disconnections as chemical reactions, in the sense that a chemical 
reaction implies an abrupt change of the objects that react. Choosing the connection as the direct reaction and 
the disconnection as the reverse reaction, we write symbolically 

v> 

(f,a) ^=» (c,a) (39) 

where the reaction rate v^ numbers the free strands which connect minus the connected strands which 
disconnect, per unit time and per unit volume in the h space. The rate v^ is assumed to obey a first order 
kinetics in the chemical sense [de Groot and Mazur (1984)], according to 

where fi c ^ and fi F ^ are the connection and disconnection rate constants, respectively, with the dimension of 
an inverse time. Since our system is assumed ideal, the rate constants depend on volume and temperature but 
not on the strand distribution, the strand being non-interacting. This excludes variation of the rate constants 
with the stress or with a structural parameter. 

In addition to the connection/disconnection reactions, free strands undergo a Brownian diffusion in the h 
space, described by the diffusion flux w£ . Its /-component w / - j . is the number of free strands having their 
connector diffusing in the /-direction of the h space, per unit area across this /-direction and per unit time. 
This Brownian diffusion is assumed to be so fast as to establish a partial equilibrium with respect to h , with 
the consequence that the chemical potential of the free strands does not depend on h , 

V F h=VF (41) 
though fi F can differ from the equilibrium value /u eq . 

The variation rate of the distribution function due to the relaxation processes can now be written 
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Fh 



dt 



~ v h 



Ch 



dt 



Wr. 
hi 



(42) 



At equilibrium, the Brownian flux and the connection rate satisfy the conditions 



eq 



vl ? =0 

h 



wi q =0 

h 



(43) 



These separate equalities, as well as the condition that vanishes (instead of the weaker conditions 

dwf jdh =0) follow from the principle of detailed balance [de Groot and Mazur (1984)]. The first line of 
eq.(43) implies that the ratio P F hJP C h cannot be chosen at will, it must satisfy the detailed balance condition 



B ~ r eq ~ = B - T eq ~ 



(44) 



The mass action law (38) allows formulating the detailed balance in terms of the connection potential: 



^Fh Fh _____ ^S/* 



J Ch 



Fh 
Ch 



■ exp 



(45) 



One more constitutive equation must be introduced, that for the heat flux q . The simplest heat conduction 
equation, Fourier's law, involves a non-negative thermal conductivity X : 



q = -XVT 



(46) 



Inserting this expression into the entropy production (16) ensures that the first term of the right member is non- 
negative. The second term can be rewritten, using eq.(40) to (42) : 



d R r 



ah 



dt 



^Fh 



Wr.- 
hi 



dhi 



+ v h A H 



(47) 



= Jrf 3 *(p cjf l 



Fh Ch 



Pi 



where Afj.~ = fj. c ~ — fi F ^ = Af- + k B T\nT c ~ jT is the chemical affinity associated to the connection reaction 
(39), and dfj.^~ jdh is the affinity of the Brownian diffusion of free strands, conjugate to the diffusion flux 
wt [Prigogine and Mazur (1953), de Groot and Mazur (1984) § X-6]. The Brownian diffusion of free strands 
gives no contribution to the source of entropy because its affinity vanishes, according to eq.(41). Using relation 
(27) and the detailed balance (45), the connection affinity can be written 
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A^=k B Tln 



Fh Ch 



(48) 



Ch Fh 



and inserting eq.(46)-(48) into the source of entropy (16) then yields 



™ = + k B \d 3 hip ci r Pi - p Pi r ci }n 



(49) 



Obviously, since the logarithm is an increasing function of its argument, the integrand is non-negative. Our 
model therefore obeys the Clausius-Duhem form (17) of the second principle : 



d irr S 
dt 



>0 



(50) 



stating that the entropy production is non-negative. This production is zero at equilibrium, when the 
temperature gradient vanishes as well as the integrated term of eq.(49), as a consequence of the detailed 
balance condition (44). Furthermore, the fact that the integrated term itself is non-negative ensures that any 
volume element d h generates a non-negative entropy [Prigogine and Mazur (1953)]. 

The distribution functions must now be related to the thermomechanical history. Using eq. (41) and (27), the 
distribution of free strands can be written as the product of the total number N F of free strands times the 
normalised distribution f ^ defined in (36) : 



N 

r Fh=J eX P 



V 



■ N ^Fh 



(51) 



The derivation of the distribution of connected strands is more involved. According to the affine deformation 
assumption, a connected strand of connector h at time t must have been connected to the network at some 
previous time f , with its connector h' related to its present value h by the affine transformation (6) : 
h' = F(t',t)-h . The probability that a strand connected at time f stays connected up to time t reads 



?(t\t',h')-- 



exp 



-JV/W(0 



with 



h"=F(t",t)h 
h' = F(t',t)h 



(52) 



where fi F ^„{t") denotes the value the disconnection rate constant had at time t" (as a function of the 
temperature and the pressure at that time). The number of free strands which connect within volume element 



d K during time interval (t',t'+dt'), and stay connected up to time t, can be written 
P' c ^,dt'V F ^, d 3 h 'p[t\t',h'), where f}' c ^, and r^, are evaluated at time ?'. Due to the affine character of the 
deformation, the volume element d h' at time t' becomes d h at time t, such that d ti = d /zdet F . 
Integrating along the past history of a given volume element d h then yields the distribution of connected 
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strands at time t : 



Y cfx = jdt'detF N' F f Fjr ,p' cjr ,p(t\t',h') 



(53) 



where h' = F h . Hereafter, all primed quantities are evaluated at the past time f. Since temperature and 
pressure vary throughout the deformation history, the temperature and pressure-dependent quantities <f>' F ^, and 
P' c ^, will vary accordingly. 

Integrating eq.(53) over h' yields an integral equation for the numbers N F and N c : 



N C =N-N F = \dt'N' F \d 3 h'(l)' F ^Pch'p(¥^') 



(54) 



Differentiating with respect to t, or equivalently integrating eq.(42) over h , gives the total connection rate: 

■ = J^Avj = N F \d 3 hP cj fi pi -\d'hH Fr T crx (55) 



dt dt 



Equations (51) and (53)-(55) permit to derive the distribution r ^ from the thermomechanical history, thus 
determining the complete set of thermodynamic functions. We now focus on the stress and the heat production. 
Inserting the distribution (53) of connected strands into expression (30) and (32) and integrating over 
d 3 h' = detFd 3 h yields 



J Chi 



dT 



hj -k B 8ij 



Jv 



(56) 



The force f c ^ = 3£ c/ - jdh and the temperature Tare of course evaluated at time t. 

Consider now the heat production. Because of its instantaneous character, the contribution of the Brownian 
diffusion of free strands to eq.(26) can be lumped into the coefficients of dT/dt and K : using eq.(51) and 
accounting for the volume and temperature dependence of <j> F ^ , the time derivative dT F ^ jdt can be written 



dV 



^Fh , dN F 
= <b ~ — — + N F 

dt Fh dt 

and the variation of entropy (18)-(22) becomes 



dS V dT V r 3 r 

— = — Cy SjjKjj - \d h 

dt T dt T v v J 



VK kk +N F 



dT 



Jv 



dT_ 
dt 



(57) 



dN r 
dt 



dT 



dT 

V J 



(58) 



where, from eq.(55), the integrated term is contributed by the connection/disconnection reactions. The quantity 
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Cw — Cy 



V 



^Fh 



dr 



Y Fh 



NAT] 



dr 



(59) 



N 



is the heat capacity at fixed distribution of the connected strands, but allowing for the fast diffusion of the free 
strands, and the tensor 



Sij=s [Y]ij +8ijN F T\d^h 



dr 



Fh 



dv 



JT 



(60) 



= ty + T8y 



dP 

jv,{T] 



■N F \d i h 



dr 



NAT) 



Y Fh 

~dv~ 



IT 



is the entropic stress tensor, also allowing for diffusion of the free strands, which contributes an isotropic term. 
The rate of heat production is obtained by inserting (58) into eq.(15) and using (41), giving 



div<? = -Cy^f + syKij +j- \d 3 h 
dt V J 



dT 



NAT) 



^ dN F * t 

dt Vfh 



dT 



NAT) 



(61) 



In deformations and temperature variations faster than the reaction rate constants f}^ , but slower than the 

Brownian diffusion of free strands, the integral term of eq.(61) becomes negligible, as well as that of eq.(58). 
Then both equations reduce to 



- . _ \dE 



(62) 



T dS 



= -C, 



dT 



V dt ' dt lJ 11 
revealing the reversible thermoelastic behaviour of the network in fast deformations. 
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5/ Entropic strands : application to polymers 



To a good approximation, polymeric strands can be modelled by entropic chains with a configurational entropy 
depending on h only, and with h -independent volume and internal energy. This situation is described by the 
chemical potential 



ah 



N 



(63) 



? ai =?„(P,r)-* Jl 7TnQ i 

where Q.^, the number of accessible configurations for a given h , is (P, 7)-independent, and £ a (P,T) does not 
depend on h . The quantity Q.^ appears in the strand entropy, given by the standard relation [de Groot and 
Mazur(1984)] 



S ah -' 



^ah 



dT 



\ dT Jp 



or 



+ k B lnH- -k B ln 



ah 



N 



(64) 



kglnD.^ being the configurational entropy. It is easy to check that expression (63) effectively leads to h - 
independent strand volume V a = (d^ /dp) = /dP^ and strand energy E a = ^ -PV a + TS aft . 
The equilibrium distribution (35) is proportional to Q.^ 

(65) 



N exp 



knT 



T Lq - = N e M r , with r = - r * and A^ q = 

J« " %' exp— ^ + exp- SF 



knT 



where 0,- is (P, 7)-independent. The detailed balance (45) then reads 



lP~ = exp 7^ ; with A Z=£c- £f 
Pen kfiT 



(66) 



and states that the ratio Pp^/Pc^ is h -independent. 

The connector force derived from expression (63) is a-independent; it derives from the configurational entropy: 

dlnQ. r 



fah~h~ k B T ~ 



dh 



(67) 



It thus satisfies the relation 



fi'T 



dT 

V Jv 



(68) 



with the consequence that that the extra stress equals the entropic extra stress 
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T ij =t ij=h\ d ^ Y cM h j- k B T5 ij) 



(69) 



In systems of en tropic strands, the free strand contribution to t is always zero because of relations (41) and 

= eq 

(68), and the whole entropic stress tensor vanishes at equilibrium, t = . 

Since the normalised distribution <j>^ is independent of pressure, volume and temperature, eq.(60) reduces to 



V 'V ~ S {T}ij ~ T H 1 



\^ JV,N, 



(70) 



i.e. the entropic stress is independent of whether the free strands are allowed to diffuse or not. A similar 
relation holds for the heat capacity (59) : 



-V{r} 



(71) 



The following relation 



dT 



dT 



(72) 



implies that \d/^ F ^ jdT ^ is h -independent like \i F ^ , therefore the Brownian diffusion of free strands 
does not contribute to eq.(22), which can be written 



— --C —- — 
dt ~ T V dt T 



V 3r JV,N a 



K kk -\d 3 hv fi 



dT 



NAT] 



(73) 



This allows rewriting the heat production equation (26) in the simple form 



div q = -C\ 



dt 



K kk +: 



dr 



-AC 



Jv 



dN r 



dt 



(74) 



When considering either a permanent network or a fast deformation, the last terms of eq.(73) and (74) can be 
neglected, and eq.(62) results in the familiar equation of rubber thermoelasticity div q =-T/V dS/dt = 

- C v dT/dt + ZyKy - T(dP/dT) v K kk [Sekimoto (1991)]. 
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6/ The limit N C «N : constant connection rate 



This section and the following one make use of the results of sections 1 to 4, but do not depend on section 5. 



The limiting case of a number of free strands N F largely exceeding the number of connected strands N c allows 
some simplifications due to the fact that N F , being close to the total number N, remains unaffected by the 
thermomechanical history of the sample. More precisely, one has N c « N = N F = N e F q and T p ^ = r e p t at 
lowest order in N c /N . The number P C ^ F f t of connection per unit time and per unit volume of the h -space, 
thus stays equal to its equilibrium value 

where g^ is called the generation function. We now rewrite the main results of section 4/ in terms of the 
generation function. The connection rate can be rewritten 

(76) 



(77) 



Y h %h $Fb^Ch ^Fh\^ch ^Ch 

The detailed balance (44) and the Boltzmann distribution (35) become 

°h rph Ch ^Fh r ^ 

This relation allows rewriting the connector force as 



dh /i /y ; 



Ch 



(78) 



Relations (77)-(78) represent how the kinetic quantities g~ and /} ^ are linked to the Boltzmann equilibrium 
distribution T eq ~ exp-f ? lk B T . Inserting the detailed balance (77) in the source of entropy (49) shows 



Ch "" r *Ch 
that the Clausius-Duhem inequality is satisfied: 



o < — = x^-forf + k B [d?h $ v A n - r r r 

The expression (53) of the distribution of connected strands becomes 

T c/; = Idt'detF g'^, p(t\t',h') 

— oo 

where h' = F h , and eq.(56) reads 



In 



Y ec i 

Ch 
T Ch 



(79) 



(80) 



OJ Chi 



dT 



hj - k B 8jj 



Jv 



(81) 
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II Isotropic and Gaussian systems, comparison with earlier transient network models 



In this section, we give the simpler form our results take under the most usual assumptions about the connector 
force law and the connection/disconnection rate, and then we compare our results with the corresponding ones 
found in the literature. 

First, consider ideal systems with isotropic C, a ^ , depending on the connector magnitude h = \h\ only. The 
connector force f ^ and its derivative with respect to T are collinear with h 



fah =hh l fah . fah 



<ah 



dh 



and 



J ah 

~dr~ 



-hh~ 



dT 



(82) 



Since the strand elasticity is not assumed to be entirely entropic, relation (68) does not hold in general. 
Under assumption (82), the extra stress and the entropic extra stress are symmetric tensors : 



t H =t ji = y\ d ^ h i h j T Ch h ~ 



dT 



N c 



k B T8 tj 



(83) 



(84) 



Jv 



Another usual assumption is that the connection/disconnection rate constant j8 r = /3 a /, is isotropic, depending 
on the magnitude of h only. This, added to the isotropy of £ a ? , makes the system completely isotropic. 
The most usual force law corresponds to Gaussian strands, which have a Hookean, a-independent force law: 



?„ft =C« +\xh 2 ; ] nh =%h ; and ^ 



ah 



y ah 



X 



3/2 



27ik B T 



exp- 



2k B T 



(85) 



For instance, the spring constant is % = 3k B T/ na for freely jointed entropic chains made with n segments of 
length a, at moderate extensions Qi«na). The connector force is the same for free and connected strands, thus 
the connection potential Af^ = A£ is h -independent. The extra stress (30) then takes the familiar form 



T.. =^\d 3 hT rT h i h j -^k B T8.. 

IJ y J Ch ' J y o IJ 

and the entropic extra stress (32) can be written 



(86) 



(87) 



It is easy to see in this particular example that t can be non-zero at equilibrium : just consider a purely 
energetic, hence T-independent spring constant %. Then the entropic extra stress reduces to a pressure, 
t = —Nk B TS /v . On the other hand, the entropic elasticity of the freely jointed chain makes the spring constant 
proportional to T, with t vanishing at equilibrium. 
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The transient network model presented here must now be compared with those belonging to the same lineage, 
namely that of Green-Tobolsky and Yamamoto. First, Green and Tobolsky (1946) considered a Gaussian 
generation function with an h -independent disconnection rate constant /5 F ^ = /5 F (using our notations). 
This model fits within our formalism, therefore it is thermodynamically consistent. Green and Tobolsky indeed 
showed that their model always leads to dissipation of mechanical work in closed-cycle deformations. 
Such is not the case of Yamamoto's model (1956-1958), which allows a free choice for P F ^, g^ and £ - . 
This contradicts our equations (77)-(78) which link the connector potential, and hence the connector force, to 
the kinetic quantities $ F z and g^ . Should this relation be enforced, Yamamoto's model would fit within ours. 
Lodge's model for Gaussian strands (1968) superposes a distribution of Green-Tobolsky networks, and shares 
the same consistency. In treating non-Gaussian strands, Lodge (1968) uses an inverse-Langevin connector force 
law, whereas he keeps a Gaussian generation function and an h -independent J3 F ^ = fi F (his eq.[6.2] and 
[6.3]). This contradicts our eq.(77)-(78) in the same way as Yamamoto does. 

The Tanaka-Edwards model (1992a-d) rests on the same basis as Yamamoto's except that it considers a finite 
number of strands, with a balance between free and connected strands corresponding to that of our section 4/. 
They allow a free choice of j8 f/ -, with Gaussian T e ^~ <x exp~xh 2 /lk B T and constant /5 C ^=/5 C . Their 
equilibrium distribution of connected strands T e l = P c ^T eq ^ m F is therefore non-Gaussian (eq.(2.23) of 

Cfa ^ F fa I 

Tanaka-Edwards (1992a), corresponding to our eq.(44)). Since on the other hand these authors also choose 
£ c ^=Xh 2 /lk B T as the potential for connected strands, relation (35) is violated and the Tanaka-Edwards 
model suffers from the same inconsistency as Yamamoto's. (in Yamamoto's and Tanaka-Edwards' notations, 
the connector force / = d(f>/dh derives from the " free energy per molecule " <j). Since these authors assume 
that the strand volume V c ^ does not depend on h , their <j) equals £ ^ - PV C ^ = £ c/ - +cst in our notations 
and their connector force / = d<j)/dh is equivalent to our f c ^ = d£ c ^ jdh ). 

The rheometric functions of transient network models have been investigated numerically where analytical 
solutions are not accessible. Takano (1974) considered Yamamoto's model with Hookean strands, a 
disconnection rate constant /5 F ? <^n and a Gaussian generation function, thus violating condition (78). 
Among the variants of the Yamamoto model examined by Fuller and Leal (1981), the non-preaveraged version 
of the Phan Thien-Tanner model introduced in their eq.(33-35) uses Gaussian strands, takes 
/3 p ^ = P (\ + on) , and uses a generation function which satisfies condition (77)-(78). This model, also 
considered by Hermann and Petruccione (1992), is thermodynamically consistent, in contrast with the other 
cases examined in Fuller and Leal (1981), where condition (77)-(78) is violated. Non-Gaussian strands, with 
h -independent /5 F ^ = /5 F , have been considered by Vrahopoulou and McHugh (1987) in a 
thermodynamically consistent model. It must be noted that the previous considerations are limited to a zero 
value of the slip parameter £ present in Fuller and Leal (1981) and Vrahopoulou and McHugh (1987), 
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corresponding to our affine deformation assumption. The thermodynamic intricacies introduced by the slip 
parameter are discussed by Larson (1983). 

A last remark concerns the incidence of thermodynamic consistency of the models used in numerical 
computation of complex flows. Since inconsistent models are able to extract energy out of nothing, numerical 
instabilities then can be due to the model itself rather that being caused by a failure of the numerical scheme. 
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Summary 



The present work applies the thermodynamic ideas of Prigogine and Mazur (1953) to the network models of 
Yamamoto (1956) and Tanaka - Edwards (1992). The formalism of chemical thermodynamics is thus applied 
to all changes experienced by the strands, namely the connection to the network or the disconnection from it, 
the variation of the connector due to the macroscopic deformation of the sample, and the Brownian diffusion of 
the free strands. The strand chemical potential relates the variation of the strand distribution to the entropy, 
energy and volume variation through the Gibbs equation. Therefore, the thermodynamic potentials keep their 
meaning when the system is driven out of equilibrium by thermal and/or mechanical action. Singling the part 
due to the deformation out of the variation of the strand distribution makes the stress naturally appear as the 
deformation thermodynamically conjugate quantity. This also introduces a new quantity we call the entropic 
stress, construed as the entropic part of the stress, relating the entropy to the deformation. The entropic stress 
appears in the heat equation, giving the heat generated or absorbed as a function of the thermomechanical 
history. The ability to cope with the thermal as well as the mechanical aspects of rheothermodynamic processes 
is an achievement of the present formalism. 

The model is kept as simple as possible : The strands form an ideal system (in the chemical sense), the 
connection/disconnection reaction obeys a first order kinetics (still in the chemical sense), the Brownian 
diffusion of free strands is fast enough for resulting in a partial equilibrium, and heat conduction is governed by 
the Fourier law. The condition that the connection/disconnection reaction stops at equilibrium introduces the 
detailed balance condition relating the reaction rate constants to the equilibrium distribution, hence to the 
connector force law. The entropy source is shown to satisfy the all-important Clausius-Duhem relation : it is 
non-negative for all possible thermomechanical histories, and it vanishes at equilibrium. This ensures the 
thermodynamic consistency of our model. 

The important special case of entropic strands leads to the expected identity between the extra stress and the 
entropic extra stress. Consideration of a constant connection rate permits to introduce the generation function 
used in most of the transient network literature. The detailed balance, the source of entropy and the stress are 
rewritten in terms of the generation function. 

A review of Yamamoto' s model progeny shows that all cases of thermodynamic inconsistency can be traced to 
the failure to obey the detailed balance condition, in the form (77) for the models expressed in terms of a 
generation function and in the form (45) for those explicitly considering the free strands. 
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Symbols, equation of definition 



a 


(3) 


A 


(20) 


Pah 


(40) 


C V 


(58) 


L ah 


(3) 


Cv,{T] 


(18) 


Y eq ~ 
ah 


(35) 


q (as subscript) (3) 


^h 


(38) 


d\d R 


(7) 


AC 


(66) 


d irr 


(15) 




(47) 


E 


(3) 




(63) 




(67) 


^ah 


(27) 


f ah' f ahi 


(29) 


e 


(27) 




(4) 


K a 


(5) 


F (as subscript) (3) 


X 


(46) 


8h 


(75) 


»ah 


(3) 


h, \ 


(3) 




(33) 


k B 


(27) 




(33) 


N 


(1) 




(13) 


N a 


(1) 




(10) 


Aff 


(36) 


= e 1 T eq 


(34) 


P 


(3) 




(65) 


p(t\t',h') 


(52) 




(36) 


q 


(14) 


X 


(85) 




(22) 


h 


(63) 


I Sy 

S ah> S » 
t, tjj 

s 

T 

v h 
V 

h hi 


(58) 

(64) 

(24) 

(3) 

(3) 

(40) 

(3) 

(42) 
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